Magnetically tunable Feshbach resonances in ultracold Li-Yb mixtures 
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We investigate the possibility of forming Li+Yb ultracold molecules by magnetoassociation in 
mixtures of ultracold atoms. We find that magnetically tunable Feshbach resonances exist, but are 
extremely narrow for even-mass ytterbium isotopes, which all have zero spin. For odd-mass Yb 
isotopes, however, there is a new mechanism due to hyperfme coupling between the electron spin 
and the Yb nuclear magnetic moment. This mechanism produces Feshbach resonances for fermionic 
Yb isotopes that can be more than 2 orders of magnitude larger than for the bosonic counterparts. 
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It has recently become possible to convert gases of 
alkali-metal atoms coherently into molecules in their ab- 
solute ground state, by magnetoassociation followed by 
stimulated Raman adiabatic passage (STIRAP). This has 
so far been achieved for KRb [1] and Cs2 [2], but other 
species are under active investigation. Other methods for 
producing ultracold ground-state molecules are also un- 
der development [5HS]. Ultracold polar molecules have 
important potential applications in quantum comput- 
ing and quantum simulation, and dipolar quantum gases 
promise a rich new physics with many novel properties 

Even more new properties would be accessible for ul- 
tracold molecules with both electric and magnetic dipole 
moments [TO]. However, alkali-metal dimers have sin- 
glet ground states, which cannot be tuned magnetically 
except through the very small magnetic moments of 
the constituent nuclei. There is therefore great interest 
in producing ultracold polar molecules with doublet or 
triplet ground states. Among the most promising candi- 
dates for this are molecules formed from an alkali-metal 
atom and a l S atom such as an alkaline earth or Yb. 
A number of experimental groups are studying ultracold 
mixtures in which such molecules might be formed, in- 
cluding Rb+Yb [H] and Li+Yb [TUH]. 

It was originally believed that systems of this type 
would not have Feshbach resonances suitable for mag- 
netoassociation, because the bound-state channels (cor- 
relating with the upper hyperfme state of the alkali- metal 
atom) are uncoupled to the entrance channel by the con- 
ventional collision Hamiltonian. However, we have re- 
cently shown that a coupling does exist, due to the mod- 
ification of the alkali-metal hyperfme coupling by the sin- 
glet atom at short range [J5]. For the prototype system 
Rb+Sr, we showed that, under favorable circumstances, 
this mechanism could produce magnetically tunable Fes- 
hbach resonances with widths up to 100 mG, which might 
be used for magnetoassociation. 

In the present paper we consider the topical system 



Li+Yb, which is under experimental investigation by at 
least two groups and has also been investigated using 
electronic structure calculations [THl HZ] • We show that 
the mechanism we previously proposed does produce Fes- 
hbach resonances in this system, but that for bosonic 
(spin-zero) isotopes of Yb they are extremely narrow. 
However, for fermionic isotopes of Yb there is a different 
mechanism that can produce resonances that are more 
than 2 orders of magnitude wider. Furthermore, the re- 
cent measurements of the scattering length for 6 Li- 174 Yb 
[T31 [2] allow us to estimate the magnetic fields at which 
these Feshbach resonances will occur. 

The collision Hamiltonian for a pair of atoms a and b 

is 
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H a + H b + V{R), (1) 



where r is the internuclear distance, /i is the reduced 
mass, L 2 is the angular momentum operator for mechan- 
ical rotation of the atoms about one another, H a and Hb 
are the Hamiltonians for the free atoms (in an applied 
field) and V(r) is the interaction operator. For collision 
of an alkali-metal atom a with a closed-shell ( 1 S) atom b, 
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where £ a is the hyperfme coupling constant for atom a, 
I, i a and ib are the electron and nuclear spin operators, 
flai 9a, uc and g" uc are the corresponding g-factors, and B 
is the magnetic field, whose direction defines the z-axis. 
There is only one electronic state, of symmetry 2 S + , aris- 
ing from interaction of atoms in 2 S and 1 S states. As de- 
scribed in ref. |15| , the hyperfme coupling constant of the 
alkali metal is modified by the presence of the closed-shell 
atom, ( a {R) = Ca + A£ a (i?). However, if the closed-shell 
atom has non-zero nuclear spin ib, there may also be a 
coupling Cb(R) — ^Cb(R) between if, and the electron 
spin; this term has not been considered previously. The 
interaction operator is therefore 

V(r) = V(R) + ACa(i?K • s + AC b (R)i b ■ s. (3) 
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FIG. 1: Interaction potential V(R) for LiYb (lower panel) 
|17| and R-dependence of hyperfine coupling constants At^(R) 
(upper panel). 



There are additional small contributions to V(R) from 
nuclear quadrupole coupling, electron-nuclear dipolar 
coupling and spin-rotation interactions, but the AC terms 
have by far the largest effect [To] , 

In the present work, we obtained the quantities CLi(i?) 
and (^Yb(R) with the relativistic density- functional the- 
ory (DFT) approach [18] implemented in the ADF pro- 
gram Ql], using the PBEO functional [U] with a QZ4P 
quadruple-zeta relativistic basis set [2"T] . The asymp- 
totic value of ^Li for 7 Li was underestimated by 5.8% 
in the DFT calculations, so we scaled (u(R) to repro- 
duce the experimental atomic value. The resulting quan- 
tities A£(i?) are shown in Figure flj they were fitted to 

— t — \ 2 — 
Gaussian forms ^oe~ a(r ~ rc > , with parameters Co = —129 

MHz, a = 0.101 A" 2 and r c = 2.60 A for 7 Li and 

Co = -406 MHz, a = 0.31 A~ 2 and r c = 3.39 A for 

173 Yb. This corresponds to a 32% maximum reduction 

in C for 7 Li (and also for 6 Li). The values for 6 Li and 

171 Yb are obtained by scaling according to the nuclear 

magnetic moments, giving Co = —48 MHz and —1500 

MHz respectively. 

Zhang et al. [T7] have carried out detailed studies of 
the potential curves of the ground and excited states of 
LiYb. In the present work, we used their CCSD(T)/ECP 
potential points [55], obtained from large-basis coupled- 
cluster calculations using a fully relativistic effective core 
potential for Yb and including counterpoise corrections. 
These were extrapolated at long range using the form 
V{r) = -C 6 R- e , with coefficient C 6 = 1594 E h a% 
[T7] , The short- and long-range parts of the potential 
were smoothly connected between 20 and 25 ao with the 
switching function used by Janssen et al. [55] , 

For a potential with known long-range behavior, both 
the energies of high-lying bound states [24] and the scat- 
tering length [25] may be expressed in terms of the frac- 
tional part of the quantum number at dissociation, or 
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FIG. 2: (Color online) Upper panel: Molecular levels (dotted 
and dashed) that cross atomic thresholds (solid) for 6 LiYb, 
calculated for potentials V+ (dotted lines) and V- (dashed 
lines) as a function of magnetic field. The Yb isotope cho- 
sen makes little difference on the scale of the Figure. The 
lower panel shows an expanded view of the crossing region 
for potential V+. Filled circles show positions where narrow 
resonances occur for all isotopes and open circles show where 
wider resonances occur for Yb isotopes with non-zero spin. 



equivalently in terms of a semiclassical phase integral, 
$ = J k(R)dR, where k 2 (R) = 2fiV(R)/h 2 . The poten- 
tial curve of Zhang et al. [T7] supports 24 bound states 
for 6 LiYb. However, it is probably accurate to only a 
few percent, and variations within this uncertainty are 
enough to span the entire range of scattering length a 
from +oo to — oo. However, Ivanov et al. [13] and Hara 
et al. [m have recently measured values of the scatter- 
ing length to be \a\ — 13 ± 3 ao (0.69 ± 0.16 nm) and 
1.0 ± 0.2 nm, respectively. Because the sign of a is un- 
known, we have scaled the potential of Zhang et al. to 
give two potentials, V+ and V-, with scattering lengths 
a = +0.8 nm and —0.8 nm, each within the error bars of 
both measurements. 

Zero-energy Feshbach resonances of the type that 
might be suitable for magnetoassociation occur at mag- 
netic fields where molecular levels cross atomic thresh- 
olds. In a molecule such as LiYb, the entire molecular 
Hamiltonian is nearly diagonal in a basis set of field- 
dressed atomic functions, so that the molecular levels lie 
parallel to the atomic thresholds as a function of mag- 
netic field. This is shown in Fig. [2j for the V + potential, 
there is a group of crossings around 1000 G that might 
in principle cause Feshbach resonances, whereas for V- 
the corresponding crossings occur around 1600 G. 

In the present work, we characterize resonances by car- 
rying out coupled-channel scattering calculations using 
the MOLSCAT program [2J], modified to handle colli- 
sions of atoms in magnetic fields [27] [28] . The molecu- 
lar wavefunctions are conveniently expanded in an un- 
coupled basis set |sm s )|iLiTOj I Li)|*Yb' 7l i,Yb)|i^L)- The 
Hamiltonian is diagonal in the total projection quan- 
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turn number M to t = Mp + Ml, where Mp = m s + 
n^i.Li + "T-i,Yt>- This provides the scattering length a(B), 
which (in the absence of inelastic scattering [25]) follows 
the functional form a(B) = Ob g [1 + A/(B — B Tes )]. The 
width A is conveniently obtained from the difference be- 
tween the positions of the pole and zero in a(B). In the 
present work, we extended MOLSCAT to provide an op- 
tion to converge on poles and zeroes of a{B), instead of 
extracting them from a fit to a grid of points. 

There are only very small terms in the collision Hamil- 
tonian that couple states of different L or Ml- For mix- 
tures containing spin-zero isotopes of Yb (mass num- 
bers 168, 170, 172, 174 and 176), conservation of M tot 
and Ml requires that the Li projection quantum num- 
ber m/^i = m a + m i,Li is also conserved. Because of 
this, Feshbach resonances can occur only where states 
of the LiYb molecule cross atomic thresholds with the 
same value of my Li, as shown by solid circles in Fig. 
[2j The lowest-field resonances for the V± and V- po- 
tentials for 6 Li 174 Yb are shown in Table [jj However, in 
both cases the resonances are extremely narrow (A < 10 
fj,G) and would be very challenging to observe experi- 
mentally. This arises mostly because the 6 Li hypcrfinc 
coupling constant is only 152 MHz in the free atom, so 
that even a 32% modification of it causes only very weak 
coupling. 

The resonance positions for 6 LiYb depend only very 
weakly on the Yb isotope, because changing the Yb iso- 
topic mass has very little effect on the reduced mass for 
collisions with a light atom such as Li. The 6 LiYb res- 
onance positions are also quite insensitive to the quality 
of the potential curve. They do depend somewhat on the 
value of Cq: for example, increasing C% by 1% shifts the 
calculated resonance positions to higher field by 0.5 G for 
fixed a. However, the main source of uncertainty in the 
6 LiYb resonance positions comes from the uncertainty in 
the measured scattering length for 6 Li 174 Yb, and hence 
on the potential scaling, as described below. 

Table [l] also includes predictions for 7 LiYb, where the 
resonances are a little wider because of the larger value 
of ( - However, in this case the 17% mass scaling in- 
creases the number of bound states from 24 to 26: the 
scattering length obtained for 7 LiYb, and the positions 
of the resonances, are far more sensitive to the depth of 
the potential. For 7 LiYb, the predicted resonance posi- 
tions could easily be in error by several hundred Gauss. 
An experimental measurement of the scattering length 
for 7 LiYb would be needed to pin this down more ac- 
curately. Nevertheless, it is clear from Table [l] that the 
Feshbach resonances are very narrow for mixtures of 7 Li 
with even-mass Yb isotopes, and this conclusion will not 
be altered by changes in the potential curve. 

There is however a different mechanism that can create 
Feshbach resonances for Yb isotopes with non-zero spin 
( 171 Yb, i = 1/2 and 173 Yb, i = 5/2). Figure [l] shows 
that, as Li approaches Yb, it polarizes the spin density of 
Yb and creates a hyperfine coupling between the molec- 
ular spin and the nuclear magnetic moment of Yb. The 



TABLE I: Calculated properties of LiYb Feshbach resonances 
caused by low-field-seeking bound states (rn^ ound — + |) 
crossing the lowest m s = —1/2 threshold state. The bound 
states are labeled by the total angular momentum projection 
of the lithium atom alone, m h f °^ d . The results for Li 174 Yb 
are for M tot = m/ jLi and those for 6 Li 173 Yb and 7 Li 173 Yb are 
for Mtot = 1 and +3/2, respectively. 



system 


-Bres (G) 


o bg (nm) 


A ( M G) 


bound 
m /,Li 


V+ potential 


6 Li 174 Yb 


996 


+0.80 


1.89 


+ 1/2 


7 Li 174 Yb 


261 


+1.74 


1.64 


+ 1 


6 Li 173 Yb 


961 


+0.82 


824 


+3/2 




986 


+0.82 


1.03 


+ 1/2 




1013 


+0.82 


<0.1 


-1/2 


7 Li 173 Yb 


213 


+1.80 


178 


+2 




253 


+1.80 


38.5 


+1 




307 


+1.80 


1.17 







388 


+1.80 


<0.1 


-1 


V- potential 


6 Li 174 Yb 


1648 


-0.80 


-6.13 


+ 1/2 


7 Li 174 Yb 


679 


+0.73 


14.9 


+ 1 


6 Li 173 Yb 


1608 


-0.76 


-2790 


+3/2 




1634 


-0.76 


-4.59 


+ 1/2 




1661 


-0.76 


< l-o-il 


-1/2 


7 Li 173 Yb 


615 


+0.76 


857 


+2 




669 


+0.76 


31.1 


+1 




733 


+0.76 


0.52 







810 


+0.76 


<0.1 


-1 



hyperfine coupling produced in this way is substantially 
larger than the change in the Li hyperfine coupling. The 
operator A^Yb(-R)*Yb • s has matrix elements that con- 
nect basis functions with Arn s = ±1 and Am^yb = +1- 
This coupling creates Feshbach resonances where states 
of the LiYb molecule cross atomic thresholds with values 
of mjLi that differ by ±1, as shown by open circles in 
Fig. | 

We have further extended MOLSCAT to handle the 
additional basis sets and couplings required for the 
ACyb(-R) coupling term, and have carried out coupled- 
channel calculations of the resonance positions and 
widths for 6 Li 173 Yb and 7 Li 173 Yb. The results are given 
in Table HJ It may be seen that there are now resonances 
more than 2 orders of magnitude wider than for Li 174 Yb, 
driven by a direct coupling involving A£yb(r)*Yb ' The 
calculated widths range from 0.8 to 2.8 mG in width for 
6 Li 173 Yb. 

The criterion for producing molecules at a narrow res- 
onance is that the field must vary smoothly enough to 
achieve adiabatic passage. Field calibration and field in- 
homogeneity are not necessarily a problem; it is high- 
frequency field noise that must be reduced to well below 
the resonance width. 
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FIG. 3: (Color online) Positions and widths of the reso- 
nances at the lowest Mtot(m s = —1/2) threshold for 6 Li 174 Yb 
(dashed) and 6 Li 173 Yb (solid), as a function of the back- 
ground scattering length produced by scaling the potential. 
The width is not shown for the very narrow m^ ^ 1 " 1 = — 3/2 
resonance. The widths are shown as ab g A because the prod- 
uct varies smoothly whereas A itself has a pole at ab g = 0. 

The main uncertainty in the resonance positions for 
6 LiYb comes from the uncertainty in the measured scat- 
tering length. Fig. |3] shows how the resonance positions 
and widths for 6 Li 1 Yb and 6 Li 173 Yb vary with the scat- 
tering length used to scale the potential. For 7 Li 173 Yb, 
as for 7 Li 174 Yb, the resonance positions depend quite 
strongly on the quality of the potential curve. However, 
the resonance widths for all isotopes depend mainly on 
the A£(r) functions, obtained here from DFT calcula- 
tions as described above, and are fairly insensitive to the 
depth of the potential. 

At the fields where resonances occur for 6 Li 173 Yb, 
m s and m^Yb are both reasonably good quantum num- 
bers. For 6 Li 173 Yb, each line in Figure [5] represents 
several levels with different values of m^yb an d thus 
M tot . These produce resonances very close together 
(within 0.2 G), but with different widths, because the 
matrix element of iyb • s between the spin functions 
(m s = +1/2, m^Yb) an d (m s = — 1/2, m^Yb + 1) is 
|[*Yb(*Yb + 1) — m i,Yb( m i.Yb + l)] 1 ^ 2 - The resonance 
widths are proportional to the square of this. An exper- 



iment is likely to use Yb trapped in a single magnetic 
sublevel and would therefore sample only one of these 
resonances; for the lowest threshold (m/.Li = +1/2), the 
different possible values for m, jb of —5/2, —3/2, —1/2, 
+1/2 and +3/2 lead to ratios in the resonance widths 
of 5:8:9:8:5, respectively. These ratios are confirmed by 
the numerical results. The resonance widths given for 
6 Li 173 Yb in Table I, for M tot = 0, are for the widest of 
the group. There is no wide (Yb-driven) resonance at this 
threshold for m^Yb = +5/2. The Yb-driven resonances 
for 171 Yb are substantially narrower than those for 173 Yb, 
partly because 171 Yb has a smaller nuclear magnetic mo- 
ment but mostly because the matrix element of iyb ■ s is 
much smaller (1 on the scale of the ratios above). 

Resonances can also occur for Li atoms in magnetically 
excited states (m/xi = —1/2 and —3/2) at the crossing 
points shown in Fig. [2j However, the molecular states 
that are produced at these resonances are quasibound at 
energies above the lowest threshold. We have character- 
ized these states by carrying out scattering calculations 
immediately below the ra/xi = —1/2 threshold and fit- 
ting them to obtain Breit-Wigner widths T [30] . For the 
V + potential, the states with m h f ™ nd = +3/2 and +1/2 
have calculated widths T/h — 715 and 2.9 Hz, corre- 
sponding to lifetimes r = Ti/T = 0.22 and 55 ms respec- 
tively, while the state with TO^™ nd = —1/2 is too nar- 
row to characterise. At decaying resonances of this type, 
the real part a(B) of the complex scattering length a(B) 
exhibits an oscillation of magnitude a ros rather than an 
actual pole [29 . We obtain a rcs — 3.4 fm, 1.3 fim and 4.0 
/im, respectively, for these three resonances. For the two 
higher-field resonances the behavior is sufficiently pole- 
like to extract values of A, which are 824 and 1.2 /iG, 
respectively. 

In conclusion, we have investigated the possibility of 
using Feshbach resonances to produce ultracold LiYb 
molecules by magnetoassociation of ultracold atoms. The 
mechanism that we previously identified for RbSr [15] ex- 
ists for LiYb, but produces only extremely narrow Fes- 
hbach resonances in this case. However, we have identi- 
fied a new mechanism, which occurs only for Yb isotopes 
with non-zero nuclear spin ( 171 Yb, 173 Yb) that can pro- 
duce Feshbach resonances that are orders of magnitude 
wider. However, even for 6 Li 173 Yb the largest calculated 
resonance width is 2.8 mG. Achieving the field stability 
to produce molecules at such a narrow resonance will be 
experimentally challenging. 
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